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Abstract 

Influence of two types of spin structures on the form of the Fermi surface (FS) 
and a photoemission intensity map is studied for t — t' — U Hubbard model. Mean 
field calculations are done for the stripe phase and for the spiral spin structure. It 
is shown, that unlike a case of electron doping, the hole-doped models are unstable 
with respect to formation of such structures. The pseudogap anisotropies are dif- 
ferent for h— and e— doping. In accordance with ARPES data for Ld2- x Sr x CuO/\ L 
the stripe phases are characterized by quasi-one-dimensional segments of FS in the 
vicinity of points M(±7r, 0) and by suppression of spectral weight in diagonal direc- 
tion k x = k y . It is shown that spiral structures display the polarisation anisotropy: 
different segments of FS correspond to electrons with different spin polarisation. 

PACS: 71.10.Fd, 74.20.Rp, 74.20.-z 

The angular resolved photoemission spectroscopy (ARPES) is effective in the study of 
electronic structure of cuprates [1, 2]. It gives the image of the Fermi surface (FS) projec- 
tion onto a C11O2 plane. Conclusions of the early works (see [1, 2], and references therein) 
were consistent with FS of the hole type centered at point Y(it, it) of two-dimensional 
Brillouin zone (BZ). Later other versions of the FS topology have been discussed. In 
particular, a presence of an electron type of the FS have been proposed with the cen- 
ter at T(0,0) in Bi 2 Sr2CaCu20s+5 (BSCCO) [3]. Revision of the problem [4, 5] partly 
confirms the original version. At the same time for La2- x Sr x CuC>2 (LSCO) it has been 
proved a transition from h-type of the FS to e-type during transition from underdoped 
(UD) to overdoped (OD) regions of the phase diagram [6, 7]. The ARPES data gave an 
evidence of the d-wave superconducting gap and of the opening of the pseudogap (PG) in 
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UD BSCCO. Recently a bilayer splitting of bands [8, 9, 10, 11] in BSCCO and the time 
reversal symmetry breaking in UD BSCCO [12] have been discovered. 

A wide use of the photoemission intensity maps in the space k x ,k y ,u poses a problem 
of extracting the FS from the ARPES data. One aspect of the problem, calculation of 
the matrix elements is discussed in [13, 14]. However the FS topology and intensities 
of various (main or shadow) segments of the FS depend also on the spin and charge 
structures of a ground state of a system. 

The goal of the present work is a model study of the influence of periodic structures on 
the FS and the photoemission intensity based on the t-t'-U Hubbard model. Unlike the 
static structures in manganites, in cuprates one deals with rather dynamical structures 
with a time scale larger than time t > 10~ 6 -^1CT 9 sec in the /z-SR experiments. So at small 
times or for the processes with energy resolution higher than h/t the local spin structures 
may be considered in quasi-static limit. Then the question arises about a correspondence 
of the structures to the ARPES data. The experimental indications of the SDW and 
CDW structures in cuprates are the incommensurate peaks in the spin susceptibility at 
q = (ir ± 5q,ir), (ir,n ± 5q) in LSCO [15], linear structures along the Cu0 2 bonds with 
period 4a (a is a lattice constant) observed in the tunnel spectra [16], periodic chess 
structure 4a x 4a around the vortex in mixed state of BSCCO [17] etc. 

Here we study the stripe and spiral structures in the simplest mean field (MF) ap- 
proach, interpret the PG, discuss examples of structures with the different types of the 
PG anisotropy and to connect some properties of MF solutions with features observed 
in some cuprates. Earlier [18] the DDW (d- density wave) structure has been proposed 
as a possible hidden order parameter (OP) in cuprates. Here a search of possible OP is 
extended to stripe and spiral structures. Stabilization of such structures originates from 
the removing of a degeneracy of states in "hot points" - the van Hove singularities (VHS) 
or in parallel segments of the FS at the so-called nesting. 

Consider a t-t'-U Hubbard model with the zero band energies 

€k = 2t(cos k x + cos k y ) + At' cos k x cos k y (1) 

Further we put t — 1, so all energies are measured in units of t. In MF approximation 
this Hamiltonian is insufficient to describe the superconducting (SC) pairing. It may 
be supplemented by an empirical pairing interaction of electrons at neighboring centers 
deduced in more accurate approaches (for example, [19]). But here we study a normal 
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state and therefore retain the original t — t' — U model. 

A periodic 2D structure is determined by two translation vectors 

Ei = (Ei x , Ely), E 2 = (E 2x , E 2y ) (2) 

and the corresponding vectors of the inverse lattice Bi,B 2 which satisfy the relations 
EiBj = 2n5ij. (The components of Ei Bi are in units of the corresponding constants 
of the original and inverse lattices). The unit cell of the structure contains n c sites with 
coordinates j = (j x ,j y ), so that each site n = n(L,j) = (ni,n 2 ) = E1L1 + E 2 L 2 + (j x , j y ) is 
determined by the integers Li, L 2 determinig the position of the unit cell and by integers 
j — (jx,jy), fixing the site inside the unit cell. 

Let k is the quasi-momentum inside the main BZ G of a given periodic structure 
unlike the momentum k from the BZ G of the original lattice. The areas of G and G are 
restricted by conditions \kBi\ < it or \k x ( y )\ < n. The order parameters of the periodic 
structures are the mean charges and spin densities on sites of the unit cell 

1 1 

r j = TT^ < >; S m = W^ K S »ML,j) > ( 3 ) 

L l 

Here fi numerates the spin components, Nl = N/n c is a number of the unit cells, n c is 

a number of sites in it. In the MF approximation the mean energy (1) is H —< T > 

+NlU J2j( r ] — Yt^Sjii) an d the wave function is determined by the occupation of the 

one-electron eigenstates xlx °f the linearized Hamiltonian 

H Un = T + N L E{2f/r/, - 2USj, j } = £ h~ k . 

j keG 

The latter is divided into independent contributions from each reduced quasi-momentum 
k. Here rj,S^j are the operators corresponding to averages (3). In the momentum 
representation the eigenstates of are expanded over a basis set of 2n c Fermi opera- 
tors 

xl = Y,4 +Bm w m *Ah (4) 



ra,<7 



where A = 1, . . . , 2n c , Bm = Biini + B 2 m 2 and m = (mi, m 2 ) is such set of integers, that 
allows the vectors k + Bm to share all phase space G. 

The matrix of eigenvectors W ma ,\ and the eigenvalues E^ x are determined by diagona- 
lization of hi in the basis set |c~ n ): 

K L k+Bm,a J 

(hj^rn^rn'yWrn'a^x = W mat \Ej, x (5) 



Here 



{h~ k )ma,m>y = & mm ' & oo> t~ k +Bm + <f(h m ' ~~ m ) [ r ~ S ^ (o" M ) CTCT '] (6) 
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with ip(j,m) = exp — j Q )Bm\; Bm = B^mi + B 2 m2- In turn, the order parameters 
(3) are expressed as 

fo,<Vl = ^E Y, Wo,v»}sMj,m -m)W^ x Ck)W m , s ,,xCk)f(E kX - fi) (7) 

Here the Puali matrices a^, a in (9) correspond to the components S^j, Tj and / is the 
Fermi function. Equations (9) determine the self-consistent MF solution. 

The intensity of the photoemission of electron with a projection k of the momentum 
on ab plane and the energy E = hu — uj is 

I(k,u) = \M(k)\ 2 A(ku)f(u)<B>Ru,k- (8) 

It is determined by the matrix element M(k), a spectral density A(ku) and the Fermi 
function /. Usual convolution is done with the Gaussian function R^k [13] with param- 
eters, which characterize a finite resolution over k and energy. The dependence of the 
matrix element M on k was studied in [13, 14]. Here for simplicity we take a constant 
value for it, since we study an effect of the structure on the spectral density A. In the 
one-electron approximation 

Mk,<*) = jfT, E \Wrn*AVmE hx -V-Uj)6 k:k+Bm (9) 

Here A = 1, . . . , 2n c and index m = (mi,m 2 ) numerates all independent Umklapp vectors 
Bm = B11712 + B 2 m2- A standard replacement of the 5- function in (9) by a function with 
finite width Q is implied. A map of I(k x ,k y ,ou = 0) allows to visualize both the main 
and shadow segments of the FS. Though the band energies are periodic functions in k 
space ( E k+BmX = E k X for any m = (mi,m 2 ) ), the intensity (8) does not posses such 
periodicity. Therefore even for the constant matrix element M in (8) various segments of 
the FS display themselves with a different intensity due to the Umklapp processes in the 
case of SDW or CDW structure. 

The intensity map allows easily to repeat the known results [20, 21, 22, 23, 24, 25, 26] 
about changing of the FS topology with increase of the doping and opening the pseudogap 
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in UD region for the homogeneous AF solutions. The energies of the lower and upper 
Hubbard bands are 

t~ kX = At' cos k x cos k y ± yj U 2 d\ + 4t 2 (cos k x + cos k y ) 2 + const (10) 

The VHS in the DOS of the lower band correspond to the "hot points" M = (=br, 0), (0, ±tt) 
A form of the FS critically depends on a sign of t'. At t' — energy E\{k) is constant 
for all k along the magnetic Brillouin zone (MBZ) boundary At t' > the energies at 
the points M are lower than at the diagonal points. This leads to the formation of the 
known hole pockets around points (n/2, 7r/2) or the electron pockets around (tc, 0), (0, re) 
in UD case of the h- or e- doped models [21]. In the same UD region the PG is open 
at x, y directions or at diagonal directions for h- or e- doped models. Such a topology 
is determined by a profile of the energies E\(k) of lower or upper Hubbard bands as a 
function of k on the MBZ boundary. But this profiles are greatly influenced also by a 
formation of the spin or charge structures in system. 

The MF calculations allow to test a stability or an instability of the homogeneous 
AF states with respect to the formation of the periodic spin or charge structures. The 
Fig.l presents the mean energy H =< H > as a function of doping 5 = |1 — n\ for the 
normal state of the hole- and electron- doped models for a set of structures. In addition 
to homogeneous paramagnet (PM) and AF solutions the next structures are considered: 

1) The stripe structure which consists of the antiphase AF stripes parallel to the y- 
axis. A fixed width 4a of the stripe (8 sites in a unit cell of the structure) is taken and 
the domain walls are chosen as centered on bonds or on the sites. Both structures have 
a similar form of the FS, but that with the bond-centered domain walls occurs to have 
some lower energy and only it is discussed below. 

2) The spiral spin structures 

< S n >= d [e x cos Qn + e y sin Qn] (11) 

with spirality vector Q = Q x = ir(r), 1) or Q = Q xy = tt(i], n), with n ~ 0.75 — 0.8. 

3) The checkered structure with antiphase AF square domain 4a x 4a. It have a higher 
energy than the above two structures, and unrealistic shape of the FS, so the results are 
not presented here. 

4) The proposed in [18] states of the orbital antiferromagnet with staggered currents 
on a plackets. In MF approximation they do not appear for t — t' — U Hubbard model 
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without additional interactions. At hole doping and U/t > 5 there exist the MF solutions 
with staggered spin currents and rotations on 7r/2 of local spins of the neighbor sites. 
But their energies are higher than those of stripe and spiral structures and corresponding 
results are not presented here. 

Most of calculations were carried out for structures with fixed commensurate period 8a, 
though the optimal size of the AF domains or spiral vector depend on doping. For spiral 
states the dependence of Q(S) has been calculated many times. Here we are interested 
only in the specific features of the FS and the pseudogap anisotropy. 

Fig. 2 shows that at the h-doping the homogeneous AF solutions are unstable with 
respect to formation of the stripe and spiral structures. Formation of these structures 
extend the doping range where the local spin retains a nonzero value (a disappearing of 
the < S^j > corresponds to matching of the energy of some structure and that of the 
paramagnet state). At the same time at e-doping the MF energies of the above structures 
appear to be higher than the energy of the homogeneous AF solution. Large stability of 
the AF state in the e-doped models corresponds to the more wide region of the AF order 
in cuprates of the n-type than in ones of the p-type. The result about the absence of the 
stripe and spiral structures in the n-type models are consistent also with a commensurate 
peak at Q ~ (n, n) in the spin susceptibility observed in NCCO, PCCO (unlike the 
incommensurate peaks in the hole doped cuprates). 

Fig. 2a presents a map of the intensity I(k,uj = 0), Eqs.(8,9), for the stripe structure 
in the model with U/t — 6, t'/t = 0.1 at n — 0.8. In case of AF MF solution the chosen 
doping is close to the "optimal" doping at which the Fermi level passes through the VHS 
in the DOS of the lower Hubbard band. The stripes (which are parallel to an y-axis) split 
the VHS. The "hot points" M x = (ir,0) and M y = (0, ir) become non-equivalent points 
and the quasi-lD asymmetric horizontal segments of the FS are formed. Fig. 2b presents 
the intensity map symmetrized over the structures with the x- or y- stripe orientation. 
The FS form differs strongly from that for the homogeneous AF state of the model with 
t' > 0. The main difference is the absence of the FS in the diagonal direction and near 
the My. This gives an evidence that the stripes open the pseudogap in these regions of 
k. The corresponding band energies E x (k) descend below the Fermi level and the work 
function A PG = « — E x > for the electron to be removed from such k is observed as the 
PG in ARPES data. The revealed anisotropy of the FS and pseudogap may change our 
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conclusions about possible symmetry of the SC order in the striped state. The retaining 
of the quasi- ID segments of the FS only suggests the ID character of conductivity along 
the stripe direction (here y-axis). With the increase of t' up to the value 0.3 there appear 
the small additional hole pockets at around kg = (n/2,7v/2) in the intensity map. 

The anisotropy of the FS and the PG suggests the revision of the possible symmetries 
of SC order coexisting with the stripes. If the structure is symmetric with respect to 
replacement x <-> y, then the d-wave SC order < cj^ct^ >~ cosA^ — cosk y is expected. 
It provides the orthogonality of the pair function to the on-site pair function cj^cj^ which 
is suppressed by the on-site repulsion U. Taking into account a new FS view (see Fib. 2a) 
one may suppose the extended s-wave SC pair function < cLc! fc . >~ cos k x + cos k y . The 
latter might be orthogonal to the on-site pair function due to the node lines k x ± k y = 
±7r. Unlike the gain of the node lines in the angular dependence of the d-SC gap, one 
can suppose a gain of the node in an analogue of a "radial" part of the pair function. 
Verification of such hypothesis requires further calculations. 

Fig. 3a presents the one-electron energies E\{k) (eigenvalues of MF problem) as func- 
tions of quasi- momentum varying along the contour 

Y(-tt, tt) - M„(0, tt) - Y(ir, tt) - M x (ir, 0) - Y(ir, -tt). 

As expected, band energies E\{k) are periodic functions with the period n/4 on the first 
horizontal segment of contour Y — M y — Y. However, the intensity map selects only the 
energy levels which are non-shadow at given k. Fig. 3b presents such map of I{k,uj) on 
the plane k, cu for k varying along the same contour. Intersections of band energies with 
the Fermi level in the vicinity of M x correspond to quasi- ID FS of Fig. 2a. 

The above pictures may be explained by an action of the spin-dependent average 
field. The main harmonic of this field is F(n) ~ cos {Q v n + y?o} with Q v = (f]7r,7i) (here 
i] = 0.75). This field shifts upstairs the zero level e(0, it) at M y pushing it from the lower 
levels 6^ at k = (±?77r, 0) near M x . The same field descends downstairs (lower the Fermi 
level) the zero level e(n, 0) at M x pushing it from the higher zero levels e(±0.257r, n) near 

My 

The ARPES data for UD LSCO [6, 7, 27] are in qualitative agreement with the proper- 
ties of the model calculations. It was shown [6, 7, 27] a presence of two different segments 
of the FS - near the points M and at diagonal directions, systematic suppression of spectral 
weight near (7r/2, ir/2) in comparison with BSCCO or OD LSCO, the straight segments of 
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the FS near (jr, 0) with the width ~ tt/2. All these observed features have been interpreted 
as an evidence of nonhomogeneous stripe structures in UD LSCO [6, 7, 27] or the combine 
action of the stripe structure and of the order- disorder [28] . Our calculations confirm this 
interpretation. Another evidence of the stripe structure is given by the observation of the 
"incommensurate" peaks in neutron scattering at k — (ir ± 5, 7r) , (ir, n ± 5) in LSCO [15]. 

Another type of the spin and angular anisotropy of the FS manifests itself in spiral spin 
structures (11) which possess the polarisational anisotropy of the FS. Different segments of 
the FS correspond to the electrons with different preferential spin polarisation. Really, the 
average field from spiral spin structure mixes one-electron states {ci.t^l+fti} provides a 
splitting of the VHS's. But unlike the stripe structures, the occupations nk, a =< cl a c ka > 
and intensities I a (kuj) of the photoemission of the electrons with fixed polarisation a 
depends on this polarisation o. Values I a (kuo) are determined by Eqs.(8,9), but without 
summing over a in right part of (9). 

Fig. 4a presents a map of the intensity I a= ^(k,uj = 0) for the spiral state with Q = 
(0.8-7T, 7r) and for the spin polarisation a =] on the axis z normal to the plane of the 
spin rotation of the spiral structure (11). Thus a photoemission of electrons with the 
up spin polarization (a =]) corresponds to the strongly anisotropic FS. For the opposite 
spin polarisation the FS coincides with the FS of Fig. 3a, reflected in the plane x = 0. 
Fig. 3b presents the FS, symmetrized over the spins and over two types of structures with 
Q = (0.87T, 7r) Q = (n, 0.87r). In the vicinity of the point M the FS have both the 
intersection with the line M — Y typical for the FS of the h-type and the same with the 
line M — T typical for the e-type of the FS. Such double intersections have been probably 
observed in ARPES data of BSCCO [1, 2]. Direct comparison with experiment would 
require to take into account the bilayer structure of BSCCO and the corresponding band 
splitting. 

Polarisational anisotropy of the FS directly reflects the existence of the spin currents 
= — in the spiral state. According to [29] it might be a reason of the time reversal 

symmetry breaking observed in the UD BSCCO in dichroism of the ARPES signal [12]. 

In the alternative hypothesis [12, 30] this effect is explained by the specific alignment of 

the circular micro currents. 

A direct observation of the polarisational anisotropy of the FS needs the selective 

over the spin polarisation measurements of the photoemission intensity. Such selective 
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measurements of the total photoemission have been achieved in the so called "spin-orbit" 
photoemission [31]. In principle, similar selectivity is possible in ARPES also. It is 
important to continue the study of the time reversal symmetry breaking, in particular, to 
clarify, whether this effect and accompanying currents have the surface or bulk character. 
Yet we cannot answer the question: if the ground state of the BSCCO might be presented 
as a set of quasi-static domains with the spiral structure and corresponding systems of 
spin currents. 

In conclusion, the MF treatment of the normal state of the t — t' — U Hubbard model 
shows that the FS topology and the PG anisotropy depends on sign of t', on the type (e- or 
h- ) of the doping and of the spin structure. For the hole-doped models the homogeneous 
AF MF state occurs to be unstable with respect to formation of the stripe phase and 
the spiral spin structure. However, in electron-doped models the lowest energy refers to 
the homogeneous AF solution. This corresponds to a wide doping range of an existence 
of the local AF order in the n-type cuprates NCCO, PCCO and to commensurate peak 
at Q = (ir, 7r) observed in neutron scattering in these cuprates. In accordance with the 
ARPES data for the LSCO, the stripe phase is characterized by the quasi-lD behavior 
of the FS in the vicinity of the "hot points" (±7r, 0), by the pseudogap opening and by 
suppression of the spectral weight in the diagonal direction and in the direction parallel 
to stripes. Such anisotropy of FS and PG does not conform with the d- symmetry of the 
SC order. For the spiral spin structure the polarisation anisotropy of FS is revealed which 
means that different segments of the FS correspond to different spin electron polarisations. 

Work is supported by the Russian Foundation of Basic Research under Grant No. 03- 
03-32141. Authors are greatfull to V.Ya.Krivnov for stimulating discussions. 
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Fig.l. 

In the top fig. the doping dependencies of the mean energy (per one site) for the hole 
doped t — t' — U model with parameters U = 4.0, t' = 0.3 for different spin structures 
are shown: for the paramagnet state (curve PM), for the homogeneous AF state (curve 
AF), for the stripe phase with a period of 8a along the x-axis (dashed line 1), for the 
spiral states with vectors Q = tt(t], 1) (curve 2) and Q = n(rj,rj) (curve 3) at rj = 0.8. 
The bottom graphics- the energies of the same structures (with the same notations) for 
the electron-doped model. For convenience the common function F(5) = U(n 2 — 1) is 
subtructed from all energies. 

Fig.2. 

a. A map of the photoemission intensity I(k,u = 0) for the antiphase AF stripes 
parallel to the y axis, having a width 4a and the bond-centered domain walls, b. The 
same, but averaged over two stripe phases of x- and y- orientation. Parameters of model 
are U = 6, t' = 0.1. 

Fig.3. 

a. The eigenenergies E\(k) of the MF problem as functions of the quasi- momentum 
running along the contour Y(— 7r, n)-M y (0, n)-Y(n, n)-M x (n, 0)-Y(n, — ir) for the stripe 
structure with period 8a along the x-axis. b. The map of the intensity I(k, u>) for k, 
running along the same contour. The map displays the same levels E x (k), but with a 
spectral weight, determined by the structure of band states. Parameters of the model are 
the same, as in Fig.2. 

Fig.4. 

a. The map of the photoemission intensity I^(k,u = 0) for the electrons with a spin 
polarisatin a =| for the spiral structure with Q = 7r(0.8, 1). b. The same but averaged 
over two spin polarisations a =t, I and over two structures with Q = (0.8tt,tt) and 
(7r,0.87r). The model parameters are U — 6, t'=0.1. 
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